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ABSTRACT 



The Takagi function r : [0, 1] — > [0, 1] is a continuous non-differentiable function constructed by 
Takagi in 1903. This paper studies the level sets L{y) = {x : r(x) = y} of the Takagi function t(x). 
It shows that for a full Lebesgue measure set of ordinates y, these level sets are finite sets, but whose 
expected number of points is infinite. Complementing this, it shows that the set of ordinates y whose 
level set has positive Hausdorff dimension is itself a set of full Hausdorff dimension 1 (but Lebesgue 
measure zero). Finally it shows that the level sets have a nontrivial Hausdorff dimension spectrum. 
The results are obtained using a notion of "local level set" introduced in a previous paper, along with a 
singular measure parameterizing such sets. 

1. Introduction 

The Takagi function r(x) is a function defined on the unit interval x G [0, 1] which was introduced by 
Takagi 11231 in 1903 as an example of a continuous nondifferentiable function. It can be defined by 



where <C x 3> is the distance from x to the nearest integer. It has appeared in a wide variety of contexts 
in probability theory, number theory and analysis, including Bernoulli convolutions lPT4l p. 195], the 
distribution of binary digit sums ( ( tf22"Tl . t24}, El), and as fractals and functions satisfying self-similar 
analogues of the Laplace equation ([25]). 

This paper considers certain properties of the graph of the Takagi function 




(1.1) 



n=0 



Q{t) := {(x,t{x)) : < x < 1} 
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Figure 1: Ordinate level set L(y) at y = 0.5 and abscissa level set L(t(x)) at x = 0.3. 

which is pictured in Figure [T] It is well known that the values of the Takagi function satisfy < 
t{x) < g . It is also known that this graph has Hausdorff dimension 1 in M. 2 , see Mauldin and Williams 
II2T1 Theorem 7]). They add the remark that they do not know whether the 1-dimensional Hausdorff 
measure of this graph is finite or infinite ( IT2T1 p. 800]). Here we study the structure of the level sets of 
this graph. We make the following definition, which contains a special convention concerning dyadic 
rationals which simplifies theorem statements. 

Definition 1.1. For < y < | the (global) level set L(y) at level y is 

L(y):={x: t(x) = y, < x < 1}. 

We make the convention that x specifies a particular binary expansion; so each dyadic rational value 
x = p- in a level set will appear twice, labeled by each of its two possible binary expansions. 

Level sets have a complicated and interesting structure, depending on the value of y. It is known 
that there are different levels y where the level set L(y) is finite, countably infinite, or uncountably 
infinite, respectively. Concerning the size of level sets, measured by Hausdorff dimension, in 1984 
Baba [6] showed that the level set £(§) has Hausdorff dimension ^, so is uncountable. The second 
author showed ( EDI ) that the Hausdorff dimension of any level set is at most 0.699 and conjectured 
Baba's example achieves the largest possible dimension. Finally, relevant to our results here, in 2008 
Buczolich [7] proved that, in the sense of Lebesgue measure, almost all level sets L(y) are finite sets. 

The object of this paper is to study properties of "generic" level sets of the Takagi function. The 
ordinate (y-axis) notion of genericity (in the Lebesgue measure sense) is to draw an ordinate y at 
random using Lebesgue measure in [0, |], and ask for properties of such level sets that hold for a full 
Lebesgue measure of such sets L(y). The abscissa (x-axis) notion of genericity is to draw a number x at 
random in [0, 1] with respect to Lebesgue measure, and then to ask for properties of level sets L(t(x)) 
that hold for a full Lebesgue measure set of x. These two sampling methods of drawing level sets are 
pictured in Figure[T]above. We focus on ordinate notions of genericity, but also treat abscissa notions in 
order to obtain our results: we observe that these two notions of genericity give quite different answers 
concerning the structure of level sets. 

A weaker notion of "generic set" is ask for properties of level sets that are generic in the Hausdorff 
dimension sense, which means to hold for a set of (full) Hausdorff dimension 1, in the appropriate 
variable. We consider this notion in both the abscissa sense and the ordinate sense. 
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1.1. Generic level sets: Results 



Our results are obtained using methods from our paper [13, which introduced the notion of "local 
level set" of the Takagi function. We summarize relevant results in Sect. [2] and [3] Local level sets 



are sets L l ° c determined locally by combinatorial operations on the binary expansion of a real number 
x; they are closed sets and each level set is partitioned into a disjoint union of local level sets. The 
structure of local level sets L l ° c is completely analyzable: they are either finite sets or Cantor sets. 
Information of the Hausdorff dimension of such sets can be deduced from properties of the binary 
expansion of x. In llT9l we introduced a certain subset Q L C [0, 1], called the deficient digit set, made 
up of the left endpoints of all local level sets; this set parameterizes all local level sets. We showed 
that Vi L is a closed set of Lebesgue measure zero. We also introduced a new singular function p,g, the 
Takagi singular function, whose points of increase are supported on Q L . This function was used to 
show that on almost all levels y there are finitely many local level sets; howe ver t here is also a dense 
set of levels in [0, 1] which have infinitely many local level sets. In Theorem 



3.2 



below we introduce 

singular measure, the Takagi singular measure, associated to the Takagi singular function, which plays 
an important role in this paper. 

Our first result concerns the cardinality of "generic" ordinate level sets in the Lebesgue measure 
sense. 

Theorem 1.2. (Ordinate generic level sets) (1) For a full Lebesgue measure set of ordinate points 
y £ [0, |] the level set L(y) is a finite set. 

(2) For a random level set L(y) with level y drawn uniformly from y E [0, |], the expected number 
of elements in L(y) is infinite. 

We prove Theorem |1.2| in Sect. [4] and Sect. [5] This result is proved using explicit calculations of 
the Takagi singular measure of various subsets of Q, L given in the fine decomposition of the deficient 
digit set S1 L made in Sect. [4] below. These calculations make use of self-similarity properties of the 
Takagi singular measure. 

Part (1) of this result was first proved in 2008 by Buczolich Q. He proves the almost everywhere 
finiteness of level sets by a method that directly studies the graph of the Takagi function. His proof 
shows the graph Q(t) = Gi{JG r (nonconstructively) partitions into an irregular 1-set Qr and a regular 
1-set Gji, and that the irregular set Qj has y-projection of Lebesgue measure and 2; -projection of full 
measure 1. Here an irregular l-set or purely unrectifiable 1-set is a set in M 2 of Hausdorff dimension 
1 that intersects every continuously differentiable curve in a set of T^-measure zero. By Besicovich's 
theorem such a set has 1-dimensional projections of measure in almost all directions, see Falconer 
iflOl Theorem 6.1.3]. A regular l-set is a set that can be covered by countably many rectifiable curves. 



In comparison our proof of Theorem 1.2 uses the Takagi singular measure, whose support lies in Q L . 
The part of the graph Q{t) that lies above L is covered by a single rectifiable curve, the flattened 
Takagi function described in lfl9l Sect. 5] (which is of bounded variation lfl9l Theorem 5.4]), so it 
should belong to the regular part Gr of Buczolich's partition. 



We contrast Theorem 1.2 with what is known about "generic" abscissa level sets in the Lebesgue 
measure sense. In lfl9l Theorem 1.4] we showed that a "generic" local level set L l ° c obtained by draw- 
ing x with the uniform distribution on [0, 1] (Lebesgue measure) is with probability one an uncountable 
set of Hausdorff dimension 0. An immediate consequence is: For a full Lebesgue measure set of ab- 
scissa points x G [0, 1] the level set L(t(x)) is uncountable. Here we advance the following conjecture 
for abscissa level sets, concerning their Hausdorff dimension. 

Conjecture 1.3. (Abscissa generic level sets) A full Lebesgue measure set of abscissa points x E [0,1] 
have level sets L(t(x)) that are uncountable and have Hausdorff dimension 0. 
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We also prove generic results in the weaker Hausdorff dimension sense. Let dim#(r) denote the 
Hausdorff dimension of a set F. Our next two results show that the set of level sets that are large in 
the sense of having positive Hausdorff dimension are themselves generic in the Hausdorff dimension 
sense. We prove results in both the abscissa generic case and the ordinate generic case. 

Theorem 1.4. (Positive Hausdorff dimension abscissa level sets) Let T^ s be set of abscissas x G [0, 1] 
such that the level set L(t(x)) has positive Hausdorff dimension, i.e. 

Tt := {x G [0, 1] : dim^(L(r(x))) > 0}. 

Then has full Hausdorff dimension, i.e. 

d[m H (Tt) = I- (1-2) 



We do not obtain information about the Lebesgue measure of T^ s . Conjecture 1.3 above is just the 
assertion that has Lebesgue measure 0. 

Theorem |1.4| is prov ed in Section [6j It is an immediate corollary of an analogous result for local 
level sets (Theorem 6.1 1, which proves Hausdorff dimension 1 even when one restricts to x E Q L . 



We define a family of Cantor sets A27., contained in V $ s P i £l L , where r > 3 is an integer parameter, 



and show that dimj^(A2 r ) — > 1 as r — > 00. Theorem 6.1 also shows: The deficient digit set £1 has 
Hausdorff dimension 1. 

Theorem 1.5. (Positive Hausdorff dimension ordinate level sets) Let T°^ d be set of ordinates y G [0, |] 
such that the Takagi function level set L(y) has positive Hausdorff dimension, i.e. 

T% d := {y : dim H (L(y)) > 0}. 

Then r°^ d has full Hausdorff dimension, i.e. 

dim^r^) = 1. (1.3) 

Here we know that the set T°^ d has Lebesgue measure 0, a result that follows from Theorem 



1.2 



Theorem |1.5| is proved in Sect. [7] We show that the Takagi function t(x) restricted to each Cantor 
set A2 r in Sect. |6]is strictly increasing and is a bi-Lipschitz map. The result follows since bi-Lipschitz 
maps preserve Hausdorff dimension. 



1.2. Dimension Spectrum of Takagi Level Sets 

Our results establish that level sets of the Takagi function exhibit a nontrivial dimension spectrum. We 
define the dimension spectrum function for Takagi level sets to be the function 

f T {a) := dim H {y : dim H (L(y)) > a}. (1.4) 

By definition this function is a nonincreasing function of a. The second author conjectured [20], and 
de Amo et al. O proved, that there are no level sets having Hausdorff dimension exceeding \, which 
gives 

f T (a) = for a > i. (1.5) 
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We trivially have / T (0) = 1 and Theorem 1.5 of this paper establishes that lim a ^ + f T (a) = 1, by 



showing that for all sufficiently large integers r, 

1 2 log r 

/rf >l~- (1-6) 
zr Zr 

By monotonicity of this function we have f T (a) > holding on some interval [0, Qo) of positive length. 
It seems reasonable to expect that f T (a) > for < a < \. Finally, we note that the assertion that 



f T (a) < 1 holds for each a > would imply the truth of Conjecture 1.3 

One can also study analogous questions for local level sets. We define the local level set dimension 
spectrum function by 

f*(a) := dim H ({x £ Q L : dim H (L l x oc ) > a}) . (1.7) 

This function is not directly comparable with f T (a) because it samples abscissa points rather than 
ordinate points. The nonexistence of level sets of Hausdorff dimension exceeding i yields 

/ T *(a) = for a> i (1.8) 



Theorem 6.1 of this paper establishes that /*(0) = lim a _ s>0 + f*(a) = 1, by showing that, for all large 



enough integers r, 

/;(^)>i- 21og ^. (i.9) 

Zr r 

Again it follows that f*{x) > on some interval [0, ct\) of positive length. 

The notion of multifractal formalism (or thermodynamic formalism) has been introduced in con- 
nection with the Holder spectrum of points of a nonsmooth function having given Holder exponent; 
see Jaffard iTTBll . |[T6l for discussion and many references. It is predicted (under suitable hypotheses) 
that the Holder spectrum dimension f(a) of a given function is a real-analytic function of a over a 
certain range, and also exhibits a convexity property over this range. The multifractal formalism fur- 
ther relates local properties of the function (local Holder exponents) to global smoothness properties 
of the function, invoving a Legendre transform of the function f{a). Our dimension functions above 
encode rather different properties of the Takagi function-size of its level sets- but in a certain sense still 
sample local properties of this function. Therefore it seems reasonable to ask whether the associated 
dimension functions f T {a) and f*(a) might have similar analytic and convexity properties. 



1.3. Discussion 

There has been much prior study of the non-differentiable structure of the Takagi function and related 
functions according to various measures. See the papers of Allaart and Kawamura (0, H), and 
references therein. We also mention the survey paper of the first author [18]. Level sets study only one 
particular aspect of the nondifferentiability of this function. 

An interesting feature of the Takagi function is that the cardinality of "generic" abscissa and 
"generic" ordinate level sets in the Lebesgue measure sense differ drastically. This difference can 
occur because sampling a point x on the abscissa favors picking level sets which are "large." It is a 
manifestation of the non-differentiable nature of the Takagi function. This difference indicates that the 
Takagi function must (in some sense) have "infinite slope" over part of its domain. In particular t(x) 
is not a function of bounded variation. 

We raise as open problems the determination of the dimension functions f T (a) and f*(a) for 
< a < 1 and of deciding whether properties predicted by the thermodynamic formalism hold. It 
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might conceivably be true that f T (a) = f*(a) holds for < a < 1; the above discussion shows it 
holds for a > i . 

Some further improvements of the results on level sets appear in Allaart [TJ, 0. In particular 
obtains results on the detailed distribution of level sets having a given finite number of elements. 

Acknowledgments. The first author thanks D. E. Knuth for raising interesting questions about the 
Takagi function (see IPT71 Problem 82, p. 103].) We thank Pieter Allaart for allowing use to include 
his substantially simplified proof of Theorem |7.1| and for bringing the work of of Buczolich [7] to our 
attention. We thank the reviewer for helpful comments. 

2. Preliminaries: Properties of the Takagi Function 
2.1. Functional Equations 

We first recall two functional equations satisfied by the Takagi function |[T9l Lemma 2.2]. 

Lemma 2.1. (Takagi functional equations) The Takagi function satisfies two functional equations, 
valid for < x < 1, the reflection equation 

t(x) = t(1-x), (2.1) 

and the dyadic self-similarity equation 

2t(|) = t(x)+x. (2.2) 

We next formulate a local self -similarity property of the graph of the Takagi function. To describe 
it we require some functions determined by the binary expansion of x. 

Definition 2.2. Let x denote a binary expansion 

°° b 

x := Y, ^ = 0.616263..-, (2-3) 

j=l 

with each bj € {0, 1}. For each j > 1 we define the following integer-valued functions. 

(1) The digit sum function Nj(x) is 

Nj(x) :=b 1 + b 2 + --- + b J . (2.4) 

We also set 

N°(x):=j-N}{x). (2.5) 

These functions count the number of l's (resp. 0's) in the first j binary digits of x. 

(2) The deficient digit function Dj(x) is given by 

Dj(x) := N?{x) - Nj(x) = j - 2Nj(x) = j - 2(&i + b 2 + • • • + bj). (2.6) 

The function Dj(x) counts the excess of binary digits b^ = over those with b^ = 1 in the first j 
digits, i.e. it is positive if there are more 0's than l's. (Note that dyadic rationals have two different 
binary expansions, and the functions Nj(x),Nj(x), Dj{x) depend on the choice of expansion.) 



6 



The local self-similarity property of the Takagi function graph is as follows (UH Lemma 2.5(2)]). 
Lemma 2.3. (Takagi self-similarity) 

Let xq = O.6162 • • • 6 n 0°° = t^t be a dyadic rational, and parameterize the interval [J^-, as 
x = xq + Tp; for < w < 1. Then there holds 

t(x) = t(xo) + i (t(io) + D n (x ) w) . (2.7) 

That is, on the dyadic interval ^r] the graph of the function t(x) is a miniature version of its 
full graph, vertically shifted by t(xq), shrunk by a factor ™, and tilted by an additive linear factor 
^D n (x )w. 

In particular, for the case of balanced dyadic rationals, which are ones with D n (xo) = (necessar- 
ily n = 2m is even), ( |2.7| ) simplifies to 

r(x) = r(x ) + T -^, (2.8) 
which comprises only a vertical shift and shrinking of the Takagi function. 



2.2. Local level sets 

The notion of local level set L l ° c is attached to the binary expansion of an abscissa point x G [0, 1]. 
We show that certain combinatorial flipping operations applied to the binary expansion of x yield new 
points x' in the same level set. The totality of points reachable from x by these combinatorial operations 
will comprise the local level set L l ° c associated to x. 
Let a binary expansion of x G [0, 1] be given: 

x:= X^ = 0.616263..., each 6, G {0,1}. (2.9) 
3=1 

The, flip operation (or complementing operation) on a single binary digit 6 is 

6:= 1-6. (2.10) 

We associate to any binary expansion x the sequence of digit positions j at which tie-values of the 
deficient digit function Dj{x) = occur, which we call balance points; note that all such j are even. 
The balance-set Z{x) associated to x is denoted 

Z(x):={c k : D Ck (x)=0}. (2.11) 

where we define Co = cq{x) = and set cq{x) < c\(x) < C2(x) < .... This sequence of tie-values may 
be finite or infinite. If it is finite, ending in c n (x), we make the convention to adjoin a final "balance 
point" c n+ i(x) = +00. We call a "block" an indexed set of digits between two consecutive balance 
points, 

B k (x) := {bj : c k {x) < j < c k +i(x)}, (2.12) 

which includes the second balance point but not the first. We define an equivalence relation on 
blocks, written B k (x) ~ B k >(x') to mean the block endpoints agree (ck{x) = cy{x') and c k +i(x) = 
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Ck>+i(x')) and either Bk(x) = By(x') or Bk{x) = By(x'), where the bar operation flips all the digits 
in the block, i.e. 

bj h> bj := 1 - bj, c k <j<c k +i. (2.13) 

Finally, we define the equivalence relation x ~ x' on two binary expansions to mean that they have 
identical balance-sets Z(x) = Z(x'), and furthermore every block Bk(x) ~ Bk(x') for k > 0. Note 
that x ~ 1 — x; this corresponds to a flipping operation being applied to every binary digit. In ||T9l we 
showed that the equivalence relation x ~ x' implies that r(x) = t(x'), so that x and x' are in the same 
level set of the Takagi function. 

Definition 2.4. The local level set L l ° c associated to x is the set of equivalent points, 

L l ° c :={x' : x'~x}. (2.14) 

We use again the convention that x and x' denote binary expansions, so that dyadic rationals require 
special treatment. 

We recall some basic properties of local level sets ( |[T9l Theorem 3. 1, Corollary 3.2]). 

1. Local level sets L l ° c are closed sets. Two local level sets either coincide or are disjoint. 

2. Each local level set L l ° c is contained in a level set: L l ° c C L{t(x)). That is, if x\ ~ X2 then 

t{xx) = t(x 2 ). 

3. Each level set L(y) partitions into local level sets 

m= U L i oc < 2 - 15 ) 

T(x)=y 

Here Q, L denotes the collection of leftmost endpoints of all local level sets. 

4. A local level set L l ° c is a finite set if the balance-set Z(x) is finite; otherwise it is an uncountable 
perfect set (Cantor set). 

2.3. Deficient digit set Vt L 

In |[T9l we studied the set of leftmost endpoints Q L of local level sets. 

Definition 2.5. The deficient digit set £l L consists of all x such that 

°° b 

n L := {x = ]T j- : Dj(x) > for all j > 1}. 

i=i 

The deficient digit set is a Cantor-type set obtained by removing a certain countable collection of 
open intervals from the unit interval, which we describe using the following definitions. 

Definition 2.6. (1) The breakpoint set B' consists of B§ = together with the collection of all bal- 
anced dyadic rationals in Q L . These are all B' = ^Ss: that have binary expansions of the form 

B' = 0.6i62...62m-i^2m 0°° for some m > 1, 



8 



that satisfy the condition 

Dj(B') > for 1 < j < 2m- 1, and D 2m {B') = 0, (2.16) 
This condition implies &2m = 1. 

(2) The small breakpoint set B is the subset of the breakpoint set £>' consisting of Bq = plus all 
dyadic rationals in B' that satisfy the extra condition that the last two binary digits 62m- 1 = 62m = 1. 

In lfT9l we used the small breakpoint set B to label the intervals removed from [0, 1] to create the 
deficient digit set Q, L . 

Definition 2.7. For each dyadic rational B = 0.b\b2---bi0l k , k > 2 in the small breakpoint set B 
(B / B$) associate the open interval Ib '■= (x(B)~ , x(B) + ) having endpoints 

x{B)- := O-Ma-fyOl^Ol) 00 
x(B) + := 0.6 1 6 2 ...6il0 fc (00)°°, 

necessarily with k > 2. We also set Jb = (4,1)- 

The following result gives properties of the deficient digit set VL L ( [ 19, Theorem 4.6].) 

Theorem 2.8. (Properties of Deficient Digit Set) 

(1) The deficient digit set Q L comprises the set of leftmost endpoints of all local level sets. It satisfies 

n L c[o,l}. 

(2) The deficient digit sum set Q L is a closed, perfect set ( Cantor set). It is given by 

n L = [0,1)\ |J I B . (2.17) 

BeB 

where the omitted open intervals Ib, for B in the small breakpoint set, have right endpoint a dyadic 
rational and left endpoint a rational number with denominator 3 • 2 k for some k > 1. 

(3) The deficient digit set Vl L has Lebesgue measure zero. 



In lfT9l Lemma 4.5] it is shown that the value of the endpoints of the removed intervals satisfies 
x+ - x~ B = t(x(B)-) - t{x{B) + ) = ( 2 . 
so that linear interpolation of a function across the a removed interval always has slope —1. 



2.4. Takagi function on deficient digit set 

In |[T9l Theorem 4.8] we proved that the Takagi function is nondecreasing on the set \£l L . Note that 

:= {^x : x e n L } = {x 6 [0, 1] : Dj(x) > for all j > 1}, 
which shows that \VL L C Vt L . 

Theorem 2.9. (I) The Takagi function is nondecreasing on the set \$l L . 

(2) The Takagi function is strictly increasing on away from a countable set of points, which 
are a subset of those rationals having binary expansions ending in O 00 or (01)°°. For each level y the 
equation y = r{x) has at most two solutions with x G \& L . Thus if x\ < x 2 < x% all belong to 
then t(xs) > t{x\). 



This result will be used in establishing the bi-Lipschitz property appearing in Theorem 7.1 
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3. Takagi Singular Measure 



A main ingredient used in this paper will be a singular measure, which is the weak derivative of the 
Takagi singular function constructed in lfl9l . We summarize the basic facts on the Takagi singular 
function, taken from |[T9l Theorems 1.5 ]. It is pictured in Figure [2]below. 

Theorem 3.1. (Takagi singular function) The function t s (x) defined by t s {x) = t(x) + xfor x G Cl L 
is a nondecreasing function on Vt L . Define its extension to all x G [0, 1] by 

t (x) := sup{r s (xi) : x\ < x with x\ G £l L }. 

Then the function t s (x) is a monotone singular function. That is, it is a nondecreasing continuous 
function having t s (0) = 0,t s (1) = 1, which has derivative zero at (Lebesgue) almost all points of 
[0, 1]. The closure of the set of points of increase of t s (x) is the deficient digit set £l L . 




Figure 2: Graph of the singular Takagi function. 

The weak derivative of the Takagi singular function is a Radon measure fis that we call the Takagi 
singular measure, and it has the following properties. 

Theorem 3.2. (Takagi singular measure) The Takagi singular function t s (x) is the definite integral of 
a nonnegative Radon measure p,s, i.e. t s (x) = Jq d/is- The measure [is has support equal to the 
deficient digit set £l L , so that 



dp, s = l- (3-19) 
o Jn L 



[ dns= f 

JO JQ 

Thus it is singular continuous with respect to Lebesgue measure. For every Borel set K in [0, 1], 

Hs{K) =meas(T S (K)), (3.20) 
where meas(-) denotes Lebesgue measure. 



Proof. By Theorem 3.1 the function t s (x) is monotone and bounded, hence of bounded (point- 
wise) variation. It follows that its distributional derivative is a finite Radon measure, call it /x^. It is 
necessarily nonnegative since r (x) is nondecreasing (cf. ||9j Theorem 1, Sect. 5.1]). 

The support of a measure fi on the real line is the closure of the set of points of increase of fi, which 



is ft by Theorem 3.1 Since meas(r2 ) = it is a singular measure, and it is singular continuous 
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because the integral t s (x) := Jq dp,s is a continuous function of x. Finally, since [0, 1] is compact, 
any Radon measure (inner regular measure) on [0, 1] is also outer regular ([ 13, Theorem E, Sect. 52]). 
Outer regularity means that for any Borel set K, ps{K) = inf k<zu{^s{U) : U an open set}. We 
know that the Takagi singular measure of an open interval U = (0, t) is 

li S (U) := f d/j s {x) = T S (t) - r 5 (0) = r s (t) = meas((r 5 (0),r 5 (t))) = meas(r s (U)) . 



o 



where meas denotes Lebesgue measure. Outer regularity of both p$ and of Lebesgue measure now 
implies that the equality holds on every Borel set K, proving (3.20 ). □ 



4. Structure of Takagi Singular Measure 



4.1 



we define a 



We now compute values of the Takagi singular measure on various subsets of Q L In 
fine partition of Q, L and in the remaining subsections we compute the singular measures of all sets in 
this partition. 

4.1. Fine partition of deficient digit set 

We partition the set L of left endpoints of local level sets into finer pieces, as follows: 

n L = ni[jn L fin (4.i) 

in which 

■= { x e n L : Dj(x) = for infinitely many j > 1}. (4.2) 
and Qf in is its complement, 

Qj in ■= {iEfi L ; Dj(x) > 1 for all sufficiently large j}. (4.3) 

The latter set can be further partitioned into subsets labelled by elements B' in the breakpoint set B' 
consisting of all dyadic rationals B' = O.&i^-fom = such that all Dj(B') > and D2 m (B') = 0. 
To each B' € B' we associate the set 

n L (B'):={x = B' + ^: x' G V}. (4.4) 

In particular for m = we have one set B' = Bq = with Q l (Bq) = \$l L ■ 
Lemma 4.1. (Fine Partition of Deficient Digit Set) The set Qf in has a partition 

B'eB' 

with union over the breakpoint set B' . Each set Q L (B') is a closed set. 

Proof. Elements x G VL L {B') have Dj(x) > for all j > 1, D 2m {x) = D 2m {B') = 0, and 
Dj(x) > for j > 2m + 1. In particular £l L (B') C ^j, n - The sets are disjoint for different B' 
because the value 2m is uniquely determined for each element of the set Q L (B'), and this determines 
the initial digits B' uniquely. Finally we see that each element x G ^fi n nas associated to it a unique 
maximal value 2m of j such that Dj (x) = 0, (j is necessarily even) and this assigns it to a particular 

n L (B'). □ 
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4.2. Singular measure mass calculations 



The Takagi singular measure dfis is not translation-invariant. However in Theorem 4.2 below we are 
able to use its self-similarity properties to compute the /i^-measure o f cer tain sets inside Sl L , namely 
the sets Sl^ and the sets Sl L (B') in the fine partition of Q L in Lemma 4.1 where B' runs over the set 
of balanced dyadic rationals that belong to Sl L . 



Theorem 4.2. For each balanced dyadic rational B' = 0.6i62---^2m = in the deficient digit set 
SI the fine partition set Q L (B') is a closed set, and its Takagi singular measure is 



»s(n L (B')) ■, 



n L (B') 



1 



22m+l ' 



(4.6) 



Proof. We already know that ns(Sl ) = 1 via Theorem 3.2 



Claim 1. The Takagi singular measure mass of \Sl L is given by 



(4.7) 



To prove the claim, we use the self-similarity relation in Lemma |2TT| l), 



2r(-x) = x + t(x), for < x < 1. 



It x E Sl L then t(x) = t l (x) so that we obtain 



2t(^x) = x + t(x) = t s (x). 



(4.8) 



Thus if x\ < X2 with both Xi E Sl L , then 



= r s (-x 2 )-T S {- Xl ) 



r{\x 2 ) + \x^-{r{ l -x 1 )+ l -x l 

^{T S (X 2 )-T S ( X1 )) +-(X 2 -X l ) 



1 f x2 

2J X1 



1 , 

lis + 2 ( x 2 - xi) . 



(4.9) 



We may rewrite this as 



J .'"2 



MS 



.'■1 



■1'2 



I'S 



< -measQaJi,^]) 



where the la st term denotes the Lebesgue measure of the interval [xi, x 2 \. 
Now by (2. 17 1 for each m > 1 we obtain a covering of Sl L using 



^cP 2m :=[0,l)\ [J I B . 



see 

|S|<2rr 



(4.10) 



(4.11) 
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in which we remove only a finite number of the "flattened" open intervals Ib corresponding to those 
B G B (the small breakpoint set) having dyadic length at most 2m. The set Vim is a closed set 
comprised of a finite number of intervals, [xj,Xj], say, having both endpoints xj, Xj G Adding up 



the relations ( |4. 10[ > over these intervals yields 

1 



J \V 2m A JVlr, 



MS 



< -meas(7 7 2m) 



(4.12) 



in which meas("P2m) denotes the Lebesgue measure of the set Vim- Next we note that the V^m form a 
nested family V2 D V4 D V% D • • • of closed sets, with 



n 1 



m=l 



Since these sets are Borel measurable and /j,s is outer regular we have 



lim / Ms = / Ms, 



1 1 J 1 1 / ns= I I's- 

2 Pirn J 5 ^ 

cf. Evans and Gariepy [9, Theorem 1, p. 2]. Now Theorem|2.8|(2), (3) together establish that 



meas("P2m) — ^ as m — > 00. 



Thus letting m — >• 00 in (4. 12 1 yields 



Ms = 7^ / MS, 



which with J QL /Us = 1 proves Claim 1. 



Claim 2. Lef B 



1 _ _k_ 

2 2rr, 



G £>' ant/ suppose that x\ = B' + /or i = 1,2, with both x\ G SI 1 '. 77ze« 



Ms 



j-i 



2 2r, 



MS + (>2 ~ X l, 



(4.13) 



Since £?' is a balanced dyadic rational, we may deduce (4. 13 1 using the formula of Lemma 2.3 
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analogous fashion to (4.9). This proves Claim 2. 

Now we complete the proof. Claim 2 yields the formula 



-,2m 



I Ms 



MS 



< me&s([x 1 ,x 2 ])- 



(4.14) 



For any B' G B' we have 



n L (B') :={x: X = B' + with x G - flH. 
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Now we may cover the set with hV^n- We apply the approximation bound (4.14) summed up over 
all the intervals in V2n, to obtain 



>2m 



< - meas(P 2 n)- 



Letting n — > oo we deduce, using meas("P2n) — > 0, that 

o2m 



MS 



US = 7T- 



This yields, since £l L (B') := B' + ^(\tt L ), that 



1 



n L (B') 



22m+l ' 



as asserted. □. 



4.3. Singular measure mass of fl^ 

The calculations of the last section yield the following result. 

Theorem 4.3. (Takagi singular measure: fine partition) Let us denote the Takagi singular measure. 
The sets Q^ in and are Borel sets, hence measurable. We have 

vs(n% n ) = i, (4.i5) 

which shows that 

/is(fi£,) = 0. (4.16) 
Consequently the image of this set under the Takagi singular function t s satisfies 

meas(r 5 (^)) = 0, (4.17) 

where meas denotes Lebesgue measure. 

Proof. Each set Q L (B') is closed, hence their disjoint union Qf in is a Borel set, hence is /im- 
measurable. The set Q. L is closed, hence fJ^ = S1 L \ Sl^ n is also a Borel set, hence is ^5-measurable, 
and 

ns(^) = vs(n L )-vs(n L fm )- 

(In fact one can easily show that the closure of Qj in is Q L .) 

Since fis(^ L ) = 1 the theorem will follow on showing ^s(^/j ra ) = 1- We have 

MOfin) = E VS(V L (B)), 
BeB' 

where £>' is the breakpoint set. Theorem |4.2| now gives /is(S7 L (£?)) = 22 „ 1 1+1 , where B = 0.b\ ■ ■ ■ &2m = 
gjlsr, with k odd. Recall from |fl9l Lemma 4.2] that the number of balanced dyadic rationals in £l L hav- 
ing the form for an odd k is the m-th Catalan number C m = ^ ( ™ ) . Here for m = we have 
Cq = 1 corresponding to the element Bq = 0. 
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The Catalan numbers are well known to have the generating function 

1 - Vl - 4z 2 



i=o 



2m 



2z 2 



(4.18) 



Inconsequence, taking z = obtain Y^jLo = ^- Therefore we obtain, using Theorem 



4.2 



that 



m=0 \m=0 



Cm 
22m 



1, 



which proves (4.15 1. Now (4.16 1 follows, and (4.17 ) follows from Theorem 3.2 on taking K = in 
d3T20]>. □ 



5. Cardinality of Global Level Sets 



In this section we prove Theorem |1.2| which states that for a full measure set of ordinates y the level 
set L(y) is a finite set, and that the expected number of elements in this set, with respect to Lebesgue 
measure on < y < | , is infinite. 

We use the following result |fl9l Theorem 5.8] giving the expected number of local level sets in a 
uniformly chosen level y in [0, |]. 

Theorem 5.1. (Expected number of local level sets) For a full Lebesgue measure set of ordinate points 
y G [0, |] the number N loc (y) of local level sets at level y is finite. Furthermore 



[ 3 N loc (y)dy = l. (5.1) 
Jo 



That is, the expected number of local level sets on a randomly drawn ordinate level y is |. 



Theorem 5.1 was proved using the coarea formula for functions of bounded variation applied to the 
flattened Takagi function t l . 



Proof of Theorem 



1.2 



(1) Let T°^ d be the set of infinite levels, i.e. 



TZ -={y- L{y) is an infinite set}. (5.2) 



To show a full measure set of ordinates have finite level sets, we show the contrapositive, that T°^ d has 
Lebesgue measure 0. We have 

?Z d c r(n£) 1J a£ c , 

in which T(f^) := {y = t(x) : x G 0^,} detects all levels that contain at least one infinite local 
level set, and 

A-ocf := {v '■ L(y) contains infinitely many different local level sets}. (5.3) 

Now [19, Theorem 7.2 (1)] shows that A^ c has Lebesgue measure 0. Thus it suffices to prove that 
r(r2^ D ) has Lebesgue measure 0. 
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Recall that t (x) = t(x) + x, for x G Vt L . Now consider t s (x) restricted to x G Q L (B) for a 
particular B € B , the breakpoint set. We write 2? = O.&1&2 • • • &2m = 2^ where is necessarily odd. 
Then x G £l L (B) if and only if 



Lemma 2.3 then gives 



We recall from Lemma 



x = B+^, with -x G -ft L . 
2 2m 2 2 



2.1 



that 2r(|(a;)) = r(s) + x if x G whence 
2 2m+1 (r(x)-r( J B)) =r(x / )+x =t s (x'), for x' €fl L . (5.4) 

Now the linear map 

y^V :=2 2m+1 (y-r(B)) 

sends the interval [t(B), t(B)+ 22 m+i ] ont o [0, 1] and it follows from the above that it sends t(Q l (B)) C 
[r(5),r(i?) + 2 2m+i ] onto the range t s (Q l ) = [0, 1]. We see from the linearity of this map that 

r(n L (B)) = [r(B),T(B) + ^k+r], 

1 



meas it' " u 



(n L (B))) 



22m+l ' 

Adding up these contributions, the summation in Theorem [43] gives that the total Lebesgue measure in 

3 J 



y G [0, 1] covered by images of these sets (counting overlaps with multiplicity) is 



2 mea S (r(0 L ( J B))) = £ /i S (n L (5)) = L 
BeB' BgB' 

(The images have some overlap, allowing their total measure to exceed the length of the interval 
[0, |].) Viewing these points x G £l L {B) as labelling left endpoints of local level sets, this says that 
a lower bound of the total number of local level set endpoints Lebesgue-integrated over < y < |, 
counted with multiplicity, is 1. Here we did not count any local level set endpoints in r(O^) := 

r (fl L \ Qf ir ^J , coming from the image of O^, since each Q L (B) is disjoint from fi^,. Theorem j^Tj 



now gives 

2 
' 3 



N loc (y)dy = 1, 

where N loc (y) counts the number of all local level set endpoints. We have already accounted for the 
full mass of this integral above, and any omitted contribution to N loc (y) coming from T(fi^) := 

r ( n L \ fi^- ) necessarily contributes an additional nonnegative amount. Thus we may conclude that 



meas(r(^)) = 0. 



as asserted. 

(2) We aim to show that the expected size of a global level set is infinite, i.e. to show that 



\L(y)\dy = +00, 
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where \L(y)\ counts the number of elements in L(y). By the discussion above we have 

2 i 

r 3 \L(y)\dy= [ \Ht(x))\i* s (x) > £ ^TT^), (5.5) 

in which 

r(S) :=|{l<j<oo: N j (B) = 0}\. 
Here each r(.B) is finite and bounded above by m if B = O.61&2 • • • &2m- We rewrite this as 



L 



2 00 T 

3 \L(y)\dy=J2^ (5-6) 

m=0 



in which 

L m := J2 2r(B) - 
BeB',\B\=2m 

Now we observe that L m , the total number of binary sequences of length 2m having AT 2ot (JB) = 0, has 
a combinatorial interpretation as counting the number the two-dimensional lattice paths of length 2m 
starting at the origin (0, 0), taking steps either (1, 1) or (1, —1), and ending at (2m, 0). Indeed these 
paths groups into collections of paths of size 2 r under the "flipping" (reflection) operation, with each 
group containing exactly one path in B'. (See the discussion and proof in Feller [12, Theorem 4, p. 90] 
and also |[T9l Lemma 4.2].) It follows that 

L m = ( 2m ) ~ 1 2*». 
m J \/Trm 



Thus the terms in the series on the right side of (5. 61 decay like fi(^U), so the series (5. 61 diverges, 
giving the result. □. 



6. Level Sets of Positive Hausdorff Dimension: Abscissa View 

We study level sets having positive Hausdorff dimension. In the paper |[T9l Sect. 3.3]) we classified 
those local level sets containing a rational number x that are of positive Hausdorff dimension: this gives 
some explicitly determinable rational ordinates y having this property. Here we show that the set T^ of 
abscissa points in f2 L that give local level sets having positive Hausdorff dimension has full Hausdorff 
dimension 1. 

Theorem 6.1. (Local level sets of positive Hausdorff dimension) Let Ajj denote the set of x G 
such that the Hausdorff dimension of L l ° c is positive, i.e. 

: = { x e n L : dim H {L l ° c ) > 0}. (6.1) 

This set has full Hausdorff dimension, i.e. 

dim H {K L H ) = 1. (6.2) 
In particular, the deficient digit set VL L has Hausdorff dimension 1. 
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Proof. It clearly suffices to prove the first assertion. For integer r > 1 let r 2r consist of all abscissas 
x G [0, 1] that satisfy: 

(i) Dj(x) > for j ^ (mod 2r). 

(ii) D 2k r(x) = for k = 1, 2, 3, ... 

These conditions are equivalent to requiring T 2r C £l L , and that all x G T 2r have the same balance-set 
Z(x) := {A; : D k (x) = 0} = 2rN. 



Claim 1. All members x G T 2r /zave /oca/ /eve/ sefs of positive Hausdorff dimension, satisfying 

log 2 1 
log(2 2r ) ~~ 2r' 



dimH(4° C )> l0g2 ' 



Claim 1 will follow from the spacing of the balance points being an arithmetic progression. This 
makes each local level set L l " c a Cantor-like set, which has a standard tree construction covered by 
2 k intervals of length 2~ 2rk , so that dim#(Lj. oc ) = Note that the particular subintervals are cho- 
sen differently at each step, so that this is generally not a self-similar construction, but the Hausdorff 
dimension lower bound proof for the Cantor set given in Falconer iTTTl Sect. 2.3] remains valid here, 
establishing Claim 1. 

Claim 1 shows that T 2r C A^, so that 

oo 

U r 2r c A L H . 

r=l 

To complete the proof it suffices to show the sets r 2r each have positive Hausdorff dimension, which 
approaches 1 as r — > oo. 



Claim 2. For all large enough r, the set T2 r has Hausdorff dimension greater than 1 — 21 ° gr . 

Claim 2 follows by observing that r 2r is itself a self-similar Cantor set in which each block of 2r 
symbols is drawn from the set 

in 

X 2r ■■={x = 7^= O.M 2 ...6 2 r : Dj(x) > for 1 < j < 2r, D 2r (x) = 0}, 

whose Hausdorff dimension is computable by the method of Falconer ifTTl Sect. 2.3]. It is well known 
that 

1 /2r N 

K^2r = C r 



r + 1 V r 
is a Catalan number. Thus we obtain 

log C r log C r 



dinitf(T 2r .) 



log 2 2r 2r log 2 ' 



However it is well known that C r = (1 + o(l)) , as the integer r — > oo. We conclude that for large 
enough r there holds 

dim//(r 2r ) > 1 . 

r 

Claim 2 now follows, and the theorem is proved. 

Note that Clai ms 1 and 2 together imply that the local level set dimension spectrum function /* (a) 



defined in Section 



1.2 



satisfies, > 1 ^ for all sufficiently large integers r. □. 
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We have C Tf s . Theorem O then gives dim H (Tj^ s ) > dim H (A^) 



Proof of Theorem 1.4 

1. yielding dim#(Pjf ) = 1. □. 

7. Level Sets of Positive Hausdorff Dimension: Ordinate View 



Our object is to prove Theorem 1.5 which asserts that the set of ordinate levels y having dim#(L(y)) > 
has Hausdorff dimension 1. We use the result on abscissas proved in the last section (Theorem |6.1| ), 
together with a property showing that the Takagi function restricted to certain small domains in [0, 1] is 
quite well behaved, i.e. it is bi-Lipschitz map. This allows the transfer of Hausdorff dimension lower 
bounds from the abscissa case treated in Sect. |6] 

7.1. Bi-Lipschitz property of Takagi function on T 2r 

The following proof of the bi-Lipschitz property incorporates a substantial simplification of the authors' 
original argument, due to P. Allaart. 

Theorem 7.1. Let Z(x) = {j : Dj{x) = 0} be the balance-set of x £ [0,1]. For r > 1 the Takagi 
function t(x) restricted to the (compact) domain 

T 2r := {x £ n L : Z(x) = 2rN}, 

is strictly increasing and is a bi-Lipschitz. map. 

Proof. We will show that if x\ < x 2 are both in T 2r then we have the bi-Lipschitz estimates 

2 2r (x 2 - x x ) > t{x 2 ) - r(xi) > ^(x 2 - xi). (7.1) 

The lower bound shows that the Takagi function restricted to T 2r is strictly increasing. 

To prove ( |7.1| ), we first reduce to the case that x\ and x 2 have binary expansions that disagree 
somewhere in the first 2r digits. If they disagree first between the 2kr + 1 and 2(k + l)r digits, then 
we can write Xi = xq + wrtn < Wi < 1, where xq = O.&1&2 • • • &2fcr is the initial block of 2kr 
digits where they agree, and w\ < w 2 . Note that w\,w 2 G A 2r by the self-similar definition of A 2r , 
and w\,w 2 disagree somewhere in their first 2r digits. Now 

x 2~Xl = {W 2 - Wt) 



and, since D 2 k r (x\) = D 2 k r (x 2 ) = 0, Lemma 2.3 gives 



t(x 2 ) - r(xi) = ^ ( r (^2) - r(«Ji)) . 

Thus it suffices to prove ( |7.1[ ) for w\ and w 2 . 

By definition all x = 0.bib 2 ... G A 2r have D m (x) > for all m 7^ 2kr, while D 2 k r (x) = for 
all k > 1. This forces b 2 kr+i = b 2 kr+2 = and &2(fe+l)r-i = ^2(fc+i)r-2 = 1. f° r ai l k > 0. We now 
suppose x\ = O.6162... and x 2 = ^.b x b' 2 ... disagree somewhere in their first 2r digits, say at the n-th 
digit, where 3 < n < 2r — 2. Now define the dyadic rationals 

xi := 0.6x62. ..6 2 r-2l0°° = 22^1, 
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x 2 :=0.6>;..4 r „ 2 10 oo = ^ T , 
which truncate the X{ at their (2r — l)-st digits, and the dyadic rationals 

x x := 0.6 1 6 2 ...6 2r _ 2 110°° = 

... 9, 

x 2 :=0.6 1 6 2 ...6 2r _ 2 110 c 



. ■ . ■ . ' , 2k 2 + 1 



-|OC 



2 2r ' 



which truncate the X{ at their 2r-th digits. We note that L> 2r _i(xj) = 1 and Z) 2r (xj) = for i = 1, 2. 
We can now write 

Xi [2r] 

Xi = Xi + for % = 1,2, 

in which xi[2r] = 0.&2r+ifr2r+2 • • • and similarly for x 2 [2r]. Since each Xi[2r] = 0.00... and each 
X{ = 0.00... we deduce |x 2 [2r] — xi[2r]| < 4, and 

1 -n 1 1 

| >X2-X! > ( X2 - Xl )-_._. 



Now x 2 — xi > so we obtain 



1 7 1 

| >x 2 -xx > -• (7.2) 



Next we estimate the difference between t(x2) and r(xi). For the lower bound, we will relate it to 
the difference between r(x2) and t(xi). Using Lemma 2.3 and £) 2r (xj) = for j = 1, 2, we obtain 

r{x i ) = r{x i )+ K £r - » for » = 1 ' 2 - 
This yields, using < r(xj[2r]) < |, 

r(x 2 ) - r(xx) > (r(x 2 ) - r(xx)) - | • ^. (7.3) 

To estimate the right side, we will relate these quantities to r(xj). Using x, = Xj + ^f, a calculation 
using Lemma 2.3 gives 

r(xi) = r(£<) + -^—x, for i = l,2. 



Thus ( |7.3[ ) becomes 

r(x 2 ) - r(xi) > (r(x 2 ) - r(xi)) - - • ^. (7.4) 

Now comes the key observation. Both xi,x 2 G 5^ L > since D m (xi) > 1 for 1 < i < 2r — 1. By 
hypothesis xi and X2 disagree at the n-th digit, with n < 2r — 2, at which 6 n = 0, b' n = 1, and we 
define the intermediate dyadic rational 

xo :=0.b' 1 b' 2 ---b' n 0^ = ^. 

We have xo € because its binary expansion is a prefix of that of x 2 . Furthermore xo < x 2 because 
xq and x 2 disagree at the (2r — l)-st digit. Thus x\ < xq < x 2 and we may now apply Theorem 2.9 [2) 
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to conclude the strict inequality t{x\) < t(x 2 ). But r(xj) (i = 1, 2) are both dyadic rationals with 
denominator at most 2 2r_1 , so we deduce that 

r(x 2 ) -r(xi) > ^zp 
Substituting this inequality into (|7.4|) yields the lower bound 



r(x 2 ) - r(xi) > | ~. (7.5) 

This proves the function t(x) is strictly increasing on T 2r , and also gives the desired bi-Lipschitz 
estimate. 

Combining (7.5 1 with the trivial upper bound estimate t(x 2 ) — t{x\) < |, and with (7.2 1 yields 

2 2r_1 (x 2 - xi) > r(x 2 ) - r(xi) > ( X2 ~ Xl )- 

This implies ( |77T] >. □ 

Remark. The Takagi function r(ac) is not a Lipschitz map on its full domain [0, 1], nor is it a Lipschitz 
function even when restricted to the domain Or. This is because it has arbitrarily steep slopes on Q L , 
as is implicit in the singular function property. 

7.2. Hausdorff dimension of T% d 



To conclude the paper we prove Theorem 1.5 
Proof of Theorem O Let 



IT* := {y : < y < ? with dim^L(y) > 0}. 



It is well known that bi-Lipschitz maps preserve Hausdorff dimension. By Theorem 7.1 the bi-Lipschitz 
property holds for the Takagi function r restricted to the compact domain T 2r . The range of this map 
is 

f 2r := {y :y = t(x), x G T 2r }, 

which therefore satisfies 

2 log r 



dim#(r 2r ) = diniH(T 2r ) > 1 



2r 



for large enough r, as shown in Claim 2 of the proof of Theorem 6.1 
Next, Claim 1 of that proof shows that every level y £ T 2r . has 

dim^(L(y)) > dim H (4° C ) > ^- 

2r 

Combining these inequalities shows that the Takagi dimension spectrum function satisfies 

/r(-)>l-^ 
J y 2r J ~ 2r 

for all sufficiently large r. Thus we have f 2r C T% d . and 

d\m H {T% d ) > 1 



r 



Letting r — > oo gives dim^(r^ ) > 1, establishing the result. □ 
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